Transformations, Lineage Comparisons, and Analysis of Down-to-Up Protomer States of Variants of the SARS-CoV-2 Prefusion Spike Protein, Including the UK Variant B.1.1.7.

Monitoring and strategic response to variants in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represent a considerable challenge in the current pandemic and for future viral outbreaks. Mutations/deletions of the virion's prefusion Spike protein may have significant impact on vaccines and therapeutics that utilize this key structural protein in their mitigation strategies. In this study, we have demonstrated how dominant energetic landscape mappings ("glue points") based on ab inito all-atom force fields coupled with phylogenetic sequence alignment information can identify key residue mutations and deletions associated with variants. We also found several examples of excellent homology of stabilizing residue glue points across the lineages of betacoronavirus Spike proteins that we have called "sequence homologous glue points." SARS-CoV-2 demonstrates the least number of stabilizing glue points associated with interchain interactions among Down-state protomers across lineages. Additionally, we computationally studied variants among the trimeric Spike protein of SARS-CoV-2 using all-atom molecular dynamics to ascertain structural and energetic changes among variants. We examined both a theoretically based triple mutant and the UK or B.1.1.7 variant. For the theoretical triple mutant, we demonstrated through alanine substitutions that three key residues could cause the transition of Down-to-Up protomer states, where the transition is characterized by the "arm" length of the receptor-binding domain (RBD) rather than the hinge angle. For the B.1.1.7 variant, we demonstrated the critical importance of mutations D614G and N501Y on the structure and binding, respectively, of the Spike protein. We note that these same two key mutations are also found in the South African B.1.351 variant. IMPORTANCE Viral variants represent a major challenge to monitoring viral outbreaks and formulating strategic health care responses. Variants represent transmitting viruses that have specific mutations and deletions associated with their genome. In the case of SARS-CoV-2 and other related viruses (betacoronaviruses), many of these mutations and deletions are associated with the Spike protein that the virus uses to infect cells. Here, we have analyzed both SARS-CoV-2 variants and related viruses, such as Middle Eastern respiratory syndrome coronavirus (MERS-CoV), in order to understand not only differences, but also key similarities between them. Understanding similarities can be as important as differences in determining key functional features of a class of viruses, such as the betacoronaviruses. We have used both phylogenetic analysis, which traces genetic similarities and differences, along with independent biophysics analysis, which adds function or behavior, in order to determine possible functional differences and hence possible transmission and infection differences among variants and lineages.

Static all-atom energetic mappings of the SARS-Cov-2 spike protein and dynamic stability analysis of "Up" versus "Down" protomer states.

The SARS-CoV-2 virion responsible for the current world-wide pandemic COVID-19 has a characteristic Spike protein (S) on its surface that embellishes both a prefusion state and fusion state. The prefusion Spike protein (S) is a large trimeric protein where each protomer may be in a so-called Up state or Down state, depending on the configuration of its receptor binding domain (RBD) within its distal, prefusion S1 domain. The Up state is believed to allow binding of the virion to ACE-2 receptors on human epithelial cells, whereas the Down state is believed to be relatively inactive or reduced in its binding behavior. We have performed detailed all-atom, dominant energy landscape mappings for noncovalent interactions (charge, partial charge, and van der Waals) of the SARS-CoV-2 Spike protein in its static prefusion state based on two recent and independent experimental structure publications. We included both interchain interactions and intrachain (domain) interactions in our mappings in order to determine any telling differences (different so-called "glue" points) between residues in the Up and Down state protomers. The S2 proximal, fusion domain demonstrated no appreciable energetic differences between Up and Down protomers, including interchain as well as each protomer's intrachain, S1-S2 interactions. However, the S1 domain interactions across neighboring protomers, which include the RBD-NTD cross chain interactions, showed significant energetic differences between Up-Down and Down-Down neighboring protomers. This included, for example, a key RBD residue ARG357 in the Up-Down interaction and a three residue sequence ALA520-PRO521-ALA522, associated with a turn structure in the RBD of the Up state protomer, acting as a stabilizing interaction with the NTD of its neighbor protomer. Additionally, our intra chain dominant energy mappings within each protomer, identified a significant "glue" point or possible "latch" for the Down state protomer between the S1 subdomain, SD1, and the RBD domain of the same protomer that was completely missing in the Up state protomer analysis. Ironically, this dominant energetic interaction in the Down state protomer involved the backbone atoms of the same three residue sequence ALA520-PRO521-ALA522 of the RBD with the amino acid R-group of GLN564 in the SD1 domain. Thus, this same three residue sequence acts as a stabilizer of the RBD in the Up conformation through its interactions with its neighboring NTD chain and a kind of latch in the Down state conformation through its interactions with its own SD1 domain. The dominant interaction energy residues identified here are also conserved across reported variations of SARS-CoV-2, as well as the closely related virions SARS-Cov and the bat corona virus RatG13. We conducted preliminary molecular dynamics simulations across 0.1 µ seconds to see if this latch provided structural stability and indeed found that a single point mutation (Q564G) resulted in the latch releasing transforming the protomer from the Down to the Up state conformation. Full trimeric Spike protein studies of the same mutation across all protomers, however, did not exhibit latch release demonstrating the critical importance of interchain interactions across the S1 domain, including RBD-NTD neighboring chain interactions. Therapies aimed at disrupting these noncovalent interactions could be a viable route for the physico-chemical mitigation of this deadly virion.